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ABSTRACT 

A simple analytical model is used to calculate the X-ray heating of the IGM for a range 
of black hole masses. This process is efficient enough to decouple the spin temperature of 
the intergalactic medium from the cosmic microwave background (CMB) temperature and 
produce a differential brightness temperature of the order of ~ 5 — 20 mK out to distances 
as large as a few co-moving Mpc, depending on the redshift, black hole mass and lifetime. 
We explore the influence of two types of black holes, those with and without ionising UV 
radiation. The results of the simple analytical model are compared to those of a full spherically 
symmetric radiative transfer code. Two simple scenarios are proposed for the formation and 
evolution of black hole mass density in the Universe. The first considers an intermediate mass 
black hole that form as an end-product of Population III stars, whereas the second considers 
super-massive black holes that form directly through the collapse of massive halos with low 
spin parameter. These scenarios are shown not to violate any of the observational constraints, 
yet produce enough X-ray photons to decouple the spin-temperature from that of the CMB. 
This is an important issue for future high redshift 21 cm observations. 

Key words: galaxies: cosmology: theory - large-scale structure of Universe - diffuse radia- 
tion - radio lines: general - quasars: general 



1 INTRODUCTION 



. One of the most startling findings made in the last few years is 
■ the discovery of super- massive black holes at redshifts ^,5.7 with 
5_h ' black hole masses of the the order of 10 9 Afo (Fan et al. 2003 and 
2006). The origin and seeds of these black holes remain uncertain. 

Currently, there are two main scenarios for creating such mas- 
sive black holes. One is as the end-product of the first metal free 
stars (Population III stars) that have formed through molecular hy- 
drogen cooling (Abel, Bryan, Norman 2000, 2002, Bromm, Coppi 
& Larson 2002, Yoshida et al. 2003). Given the low cooling rate 
provided by molecular hydrogen, the collapsing initial cloud is ex- 
pected not to be able to fragment into small masses and thus pro- 
duce very massive stars (for reviews, see Bromm & Larson 2004, 
Ciardi and Ferrara 2005). These stars are expected to burn their 
fuel very quickly and to produce black holes with masses in the 
range 30 - 1000A/ o (O'Shea and Norman 2006), with the excep- 
tion of the mass range of 140 — 260M© where the pair-instability 
supernovae leave no black hole remnants (Bond, Arnett & Carr 
1984; Heger & Woosley 2002; Rakavy, Shaviv & Zinamon 1967). 
Such objects grew their masses very efficiently by accretion up to 
1O 9 M by z w 6 (Volonteri & Rees 2005, Rhook & Haehnelt 
2006). 

The second avenue for producing even more massive black 
holes is through the collapse of very low angular momentum gas 
in rare dark-matter halos with virial temperatures above 10 4 K (see 



Shapiro 2004 for a recent review). Under such conditions, atomic 
cooling becomes efficient and black holes with masses 3> 10 Mq 
can be formed (Bromm & Loeb 2003). Fragmentation of the initial 
gas into smaller mass objects due to efficient cooling can be pre- 
vented by trapping the Lyman-a photons within the collapsing gas 
(Spaans & Silk 2006). 

Notwithstanding the origin of these massive black holes, their 
impact on the intergalactic medium is expected to be dramatic in 
at least two ways. Firstly, these objects produce very intense ion- 
ising radiation with power-law behaviour that creates a different 
ionization pattern around them from that associated with thermal 
(i.e., stellar) sources. The ionization aspect of the miniquasar ra- 
diation has been explored by several authors (Madau, Meiksen & 
Rees 1997, Ricotti & Ostriker 2004a, 2004b, Madau et al. 2004, 
Zaroubi & Silk 2005). Recently, however, it has been argued (Di- 
jkstra, Haiman & Loeb 2004, Salvaterra, Haardt & Ferrara 2005) 
that miniquasars can not reionize the Universe as they will produce 
far more soft X-ray background radiation than currently observed 
(Moretti et al. 2003; Soltan 2003) and at the same time satisfy the 
WMAP 3 rd year polarisation results (Page et al. 2006 ; Spergel et 
al. 2006) and the reionization constraints from the IGM temper- 
ature at redshift « 3 (Theuns et al. 2002a, 2002b, Schaye et al. 
2000). It should be noted, however, that the Dijkstra et al. 2004 & 
Salvaterra et al. 2005 calculations have been carried out assum- 
ing specific black-hole mass density evolution histories and spec- 
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tral energy distributions of UV/X-ray radiation emanating from the 
miniquasars. 

Secondly, due to their x-ray radiation, even the intermedi- 
ate mass black holes (IMBH) are very efficient in heating up 
their surroundings. Nusser (2005) has pointed out that this heat- 
ing facilitates observation of the redshifted 21 cm radiation in ei- 
ther emission or absorption by the neutral hydrogen in the high- 
redshift IGM. The observation of this radiation is controlled by 
the 21 cm spin temperature, T sp i„, defined through the equation 
n\ /no = 3 exp(— r*/T 3p i n ). Here ni and no are the number den- 
sities of electrons in the triplet and singlet states of the hyperfine 
levels, and T* = 0.0681 K is the temperature corresponding to 
the 21 cm wavelength. For the 21 cm radiation to be observed rel- 
ative to the CMB background, it has to attain a different tempera- 
ture and therefore must be decoupled from the CMB (Wouthuysen 
1952; Field 1958, 1959; Hogan & Rees 1979). The decoupling is 
achieved through either Lyman— a radiation or collisional excita- 
tions and heating. For the objects we are concerned with in this 
paper, i.e., miniquasars, the collisional excitation and heating are 
much more important. In general, throughout this paper, we will 
ignore the influence of Lyman-a photons emitted by the quasar 
on T S pi n . However, one should point out that collisional excita- 
tions due to x-ray photons results in a "secondary" Lyman— alpha 
pumping which will dominate the spin temperature and CMB tem- 
perature decoupling in some regions around the miniquasar; this 
effect has been recently point out by Chuzhoy, Alvarzez & Shapiro 
(2006). For recent papers that discuss X-ray heating, see Chen & 
Miralda-Escude (2006) and Pritchard & Furlanetto (2006). 

Collisional decoupling of T 3p i n from Tc m b is caused by very 
energetic electrons released by the effect of the x-ray miniquasar 
radiation on the IGM. Shull & van Steenberg (1985) have esti- 
mated that more than a tenth of the energy of the incident photons 
is absorbed by the surrounding medium as heating (this fraction in- 
creases rapidly with the ionized fraction). The increase in the tem- 
perature is observable at radio frequencies in terms of the differen- 
tial brightness temperature, 8Tb, which measures the 21 cm inten- 
sity relative to the CMB. A similar fraction of the absorbed energy 
also goes into collisional excitation, where this fraction decreases 
rapidly with the ionised fraction. These two processes, heating and 
excitation decouple the spin temperature from the CMB tempera- 
ture and render the IGM observable through its 21 cm emission. 

Recently, Kuhlen & Madau (2005) and Kuhlen, Madau & 
Montgomery (2006) have performed a detailed numerical study of 
the influence of 150A/q IMBH on its surroundings and calculated 
the gas, spin and brightness temperatures. They have shown that 
heating by 150M© IMBH at z — 17.5 can enhance the 21 cm 
emission from the warm neutral IGM. The filaments enhance the 
signal even further and may make the IGM visible in future radio 
experiments (e.g., the LOFAR-Epoch of Reionization key science 
projecf]]). 

In this paper, we adopt a complementary theoretical approach 
to the numerical one adopted by Kuhlen & Madau (2005). This al- 
lows us to explore the influence of power-law radiation fields from 
a range of black hole masses that are presumed to reside in the 
centres of primordial miniquasars. Furthermore, the effect of X-ray 
induced collisional excitations on the 21 cm spin temperature is in- 
cluded (Chuzhoy, Alvarez & Shapiro 2006) - this effect is not taken 
into account in the Kuhlen & Madau (2005) work. We test two main 
classes of x-ray emitting miniquasars, those with UV ionizing ra- 
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diation and those without. We show that in both cases these mini- 
quasars might play an important role in heating the IGM without 
necessarily ionizing it completely. 

In addition, two simple scenarios for the formation of 
(mini)quasars as a function of redshift are presented. This is done 
using the extended Press-Schechter algorithm to predict the num- 
ber density of forming black holes either with H2 cooling or with 
atomic cooling. We also discuss the implications of these scenarios 
for the mass density of quasars at redshift 6, the soft X-ray back- 
ground (SXRB)in the energy range 0.5 — 2 keV (Dijkstra, Haiman 
& Loeb 2004), the number of ionizing photons per baryon and, fi- 
nally, the optical depth for Thomson scattering of CMB photons. 

The paper is organised as follows: Section[2]describes the the- 
oretical methods used here and derives the ionization and kinetic 
temperature profiles around miniquasars without UV ionizing ra- 
diation. Section [3] calculates the spin and brightness temperature 
around the same quasars. In section [4] we show the ionization and 
heating profiles around quasars with UV ionizing radiation. Sec- 
tion [5] present the two formation scenarios and their implications. 
The paper concludes with a summary section (§|6j. 



2 HEATING AND KINETIC TEMPERATURE 

The exact shape of the UV and X-ray photon spectral energy dis- 
tribution around high redshift mini-quasars is uncertain. However, 
in general it is believed to have two continuum components. The 
first is through to emanate from the putative accretion disk around 
a black hole; this component, at least in low mass black holes, 
is well described by "multicolour disk blackbody"(Mitsuda et al. 
1984). The hottest blackbody temperature, T ma x, in a Keplerian 
disk damping material onto a black hole at the Eddington rate 
is kT max » fclkeV(M/M )- 1/4 (Shakura & Sunyaev 1973), 
where the hole mass, M, is measured in solar mass units. The 
characteristic multicolour disk spectrum follows a power law with 
L_b1 oc E 1 ' 3 at E < kTm ax . The second component spectrum is 
a simple power law with spectral energy distribution proportional 
to E~ a with aapproxl. The precise origin of this power law is 
uncertain and very likely to be due to nonthermal processes. 

To simplify the calculation, we follow Kuhlen et al. 2005 and 
consider miniquasars with power-law flux spectra and power-law 
index of —1. We also assume, at this stage, that the ionizing UV 
photons produced by the miniquasars are absorbed by the immedi- 
ate black hole environment. Therefore a lower cutoff of the photon 
energies is assumed, namely, 

F{E) = AE^ 1 s _1 {200 eV < E < 100 keV}. (1) 

where A is normalized such that the miniquasar luminosity is a 
tenth of the Eddington luminosity. Miniquasars with UV ionizing 
photons are considered in a later stage in this paper. 

This spectrum translates to a number of photons per unit time 
per unit area at distance r from the source, 

M(E;r) = e-C^-^-iTWV 1 , (2) 
with 

r{E;r)= / n H x H icr(E)dr . (3) 
Jo 

Here xhi is the hydrogen neutral fraction, tih ~ 1.9 x 
10 -7 cm _3 (l + z) 3 (Spergel et al. 2006) is the mean number den- 
sity of hydrogen at a given redshift, and an(E) = ao (Eo/E) 3 
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is the bound-free absorption cross-section for hydrogen with ao = 
6 x 10 -18 cm 2 and Eq = 13.6 eV. The second equation is ob- 
tained assuming a homogeneous density for the IGM. 

The cross-section quoted earlier does not take into account 
the presence of helium. In order to include the effect of helium, we 
follow Silk et al. (1972) who modified the cross section to become 



a(E)=a H (E) + ^CTHe 



cri 



(§)■ 



(4) 



A proper treatment of the effect of helium is accounted for by defin- 
ing o i to be a step function at the two helium ionisation energies 
corresponding to He I and He II. This however includes lengthy cal- 
culations and complicates the treatment, and we therefore choose 
cri to be a smooth function of E, an approximation that will over- 
estimate u(E) for low energy photons. For the kinds of spectra and 
energies we consider here, this is a reasonable assumption. 

2.1 Ionization 

To obtain the optical depth at a given distance, r, from the mini- 
quasar, we calculate the neutral fraction around the miniquasar 
for a given spectrum and energy range by solving the ionization- 
recombination equilibrium equation (Zaroubi & Silk 2005): 



AV 



n H (l~x H i) =r(r) n H XHi (l + 



tTffe TflHe 



(5) 



Here V(r) is the ionisation rate per hydrogen atom at distance r 
from the source. Since we are interested in the detailed structure of 
the ionisation front, T is calculated separately for each value of r 
using the expression, 



r(r) = I a{E)N(E-r) x ( 



E , , .\ dE 



(6) 



The function 4>{E, x e ) is the fraction of the initial photon energy 
that is used for secondary ionizations by the ejected electrons and 
x e is the fraction of ionized hydrogen (Shull & van Steenberg 1985; 
Dijkstra, Haiman & Loeb 2004). The <f>(E,x e ) term is intro- 
duced to account for the number of ionization introduced by sec- 
ondary ionization. Furthermore, in eq.[5]a^ is the recombination 
cross-section to the second excited atomic level and has the values 



of 2.6 x 10 



-13/ 



5 cm s , with T4 being the gas temperature 



in units of 10 4 K. For this calculation we assume that T = 10 4 K. 

This is of course not very accurate, although it gives a lower limit 

(2) 

on the recombination cross-section, a HI (in neutral regions atomic 
cooling prevents the gas from having a higher temperature). Since 
the region we are going to explore is mostly neutral, an accurate 
estimation of the recombination cross-section is not necessary. 

Figure[TJshows the solution of equation[5]for miniquasars with 
masses ranging from 50Mq up to 2.5 x 1O 4 M0. We assume that 
the miniquasars emit at a tenth of the Eddington luminosity and 
that their emitted radiation is confined to 200 < E < 10 5 eV. 
The lack of ionising UV photons results in a very small ionized re- 
gion around the miniquasar centres (x-ray photons are not very ef- 
ficient in ionization) with an extended transitional region between 
the ionised and the neutral IGM (Zaroubi & Silk 2005). We also 
assume that the density of the IGM around the miniquasars is the 
mean density in the Universe (this could be easily replaced by any 
spherical density profile). Due to the increase of the mass density 
at higher redshifts, the ionising photons are absorbed closer to the 
quasar. The neutral fraction profile obtained for each profile is used 
in the following sections to calculate the kinetic, spin and bright- 
ness temperatures of the IGM surrounding the miniquasars. 




0.001 



10 100 1000 10 100 1000 
R [comoving kpc] R [oomoving kpc] 

Figure 1. The neutral hydrogen fraction as a function of distance for a range 
of black hole masses for z=17.5 (left) and z=10 (right) for miniquasars 
without ionising UV radiation, namely, with radiation that spans the energy 
range of 200eV < E < 10 5 eV 



2.2 Heating 

The heating rate per unit volume per unit time that is produced 
by the photons absorbed by the IGM for a given photon energy at 
distance r from the source is Ti(r). TC is calculated separately for 
each r using the expression, 



H(r) = fn H x H i(r) 



<j{E)N{E;r)dE 



(7) 



B 



where / is the fraction of the absorbed photon energy that goes 
into heating through collisional excitations of the surrounding ma- 
terial (Shull & van Steenberg 1985). The function / is fitted in 
the Shull and van Steenberg (1985) paper with the following sim- 
ple fitting formula: / = C [l - (1 - x a ) b ], where C = 0.9771, 
a = 0.2663, b = 1.3163 and x = 1 — xhi is the ionized fraction. 
This fitting function is valid in the limit of high photon energies, 
an appropriate assumption for the case at hand. We only modify 
the fitting formula by imposing a lower limit of 11% for the frac- 
tion of energy that goes into heating as the proposed fitting for- 
mula does not work well at ionized hydrogen fractions smaller than 
10 -4 . This equation is similar to that obtained by Madau, Meiksin 
&Rees(1997). 

In order to determine the temperature of the IGM due to this 
heating, we adopt the following equation, 



= H(r) x t q 



(8) 



3 n H k b T kin (r) 
2 fx 

Here, Thin is the gas temperature due to heating by collisional 
processes, ki is the Boltzmann constant, fi is the mean molecu- 
lar weight and t q is the miniquasar lifetime. This equation assumes 
that the heating rate due to the absorption of x-ray photons during 
the miniquasar lifetime is constant. Given the miniquasar lifetime 
relative to the age of the Universe at the redshifts we are inter- 
ested in, cooling due to the expansion of the Universe can be safely 
neglected. Notice that in the highly ionized regions, although we 
ignore it, Compton cooling off CMB photons for long living mini- 
quasars and high redshifts should be included. 

Figure [2] shows the kinetic temperature as a function of radius 
for the same black hole masses considered in figure[TJ The heating 
of the IGM is clearly very extended and ranges from about a quar- 
ter of a comoving Mpc for a black hole with 50 M© up to more 
than 3 comoving Mpc for black holes with masses > 10 4 M @ . 
Since the mass density in the Universe increases towards higher 2 
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Figure 2. The kinetic temperature of the gas for a range of black hole 
masses. The redshift and quasar lifetime (tq) are specified on each panel. 
The dashed line indicates the CMB temperature at the corresponding red- 
shift. 

as (1 + z) 3 , the neutral fraction around the miniquasar is larger, 
hence, the heating is more effective at higher redshifts. The figure 
also shows that, as expected, the heating is larger for a quasar with 
a longer lifetime. Note, that at redshift 10 other effects (e.g, Lyman- 
a pumping, metal cooling lines) might play a more important role 
than heating by miniquasars. However, the purpose of presenting 
the z = 10 figures is to show the redshift trend of change due to 
miniquasars. 

2.3 Comparison with a spherically symmetric full radiative 
transfer code 

In order to test our analytical approach we compare our results with 
those obtained by running a non-equilibrium spherically symmetric 
radiative transfer code that is applied to the same problem. Details 
of the code are described by Thomas & Zaroubi (2006) but here we 
give a brief description. The radiative transfer code evolves non- 
equilibrium equations for H I , H II , He I , He II , He III , e and 
the electron temperature T e . The equations take into account col- 
lisional and photo-ionization, recombination, collisional excitation 
cooling, recombination cooling, free-free cooling, Hubble cooling, 
Compton heating and Compton cooling. The comparison between 
the analytical and the numerical results is performed for 8 cases. 
The 8 cases constitute all combinations of two black hole masses 
(100 and 10000), two redshifts (10 and 17.5) and two miniquasar 
lifetimes (3 and 20 Mega-years). The comparison is shown in fig- 
ure [3] where the kinetic temperature of the gas obtained from the 
simple analytical calculation is represented by the solid line and 
that obtained from the radiative transfer code is represented by the 
dashed line. Except at the centre where the neutral fraction adopted 
profile differs in the two cases, the agreement between the two ap- 
proaches is very good. The main reason for the departure in the 



Figure 3. This figure shows a comparison between the model adopted in 
this study and the results from a spherically symmetric radiative transfer 
code (Thomas & Zaroubi 2006) applied to two of the IMBH masses, 100 
& 10000 Mq with the same radiation power spectrum. The analytical cal- 
culation is represented by the solid line and that obtained from the radiative 
transfer code is represented by the dashed line. The dotted line indicates the 
CMB temperature at the corresponding redshift. 



centre is that the equilibrium solution assumes that the neutral frac- 
tion profile shown in fig [T] is attained within the quasar lifetime; 
this assumption is simply incorrect for high energy photons where 
the bound-free time scales exceeds that (for a recipe to mitigate this 
effect see Thomas & Zaroubi 2006). To summarise, given the many 
processes included in the radiative transfer code, this agreement is 
satisfactory. 

Another comparison one can make is with the gas temper- 
atures obtained by Kuhlen & Madau (2005) shown in the upper 
right panel of figure 7 in their paper. Visual inspection of the results 
of our approach when applied to a 150 Mq IMBH with the same 
spectrum shows good agreement. Both of these comparisons give 
us confidence in the validity of the simplistic theoretical approach 
adopted in this paper. 



3 21-CM SPIN AND BRIGHTNESS TEMPERATURES 
3.1 The Spin Temperature 

In his seminal paper, Field (1958; see also Kuhlen, Madau & Mont- 
gomery 2006) used the quasi-static approximation to calculate the 
spin temperature, T ap i„, as a weighted average of the CMB tem- 
perature, the gas kinetic temperature and the "light" temperature 
related to the existence of ambient Lyman-a photons (Wouthuysen 
1952; Field 1958). The spin temperature is given by: 

rp _ T» + TcMB + VkinTkin + VoTkin 



© 0000 RAS, MNRAS 000, 000-000 



Heating of the Intergalactic Medium by Primordial Miniquasars 5 



where Tomb is the CMB temperature and ykm and y a are the 
kinetic and Lyman-a coupling terms, respectively. 

The Kinetic coupling term is due to the increase in the kinetic 
temperature due to X-ray heating. 



Vkiri 



T, 



Au)Thi 



■ (C H + C e + C p ) 



(10) 

(Wild 1952) is the Einstein 



Here A w = 2.85 x 10" lb s 
spontaneous emission rate coefficient. Ch, C e and C p are the 
de-excitation rates due to neutral hydrogen, electrons and pro- 
tons, respectively. These rates have been calculated by several au- 
thors (Field 1958; Smith 1966; Allison & Dalgarno 1969; Zygel- 
man 2005). In this paper we use the fitting formulae used in 
Kuhlen, Madau & Montgomery (2006) which we repeat here for 
completeness, the rate due to neutral hydrogen Ch = 3.1 x 
10~ 1:L 7ihT^ 357 exp(-32/r fc4n )[s" 1 ]; the rate due to electrons is 
C e — n e 7 e where log(7 e /l cm 3 s" 1 ) = —9.607 + 0.5 log T kin x 
exp (— (log Xfcj n ) 4,5 /1800) ; and the rate due to protons is C v = 
3.2 n p k, where « = Ch/tir is the effective single-atom rate 
coefficient. And tih, n e and n p are the hydrogen, electron and pro- 
ton number densities in the unit of cm -3 , respectively and T k i n is 
measured in K. 

The Lyman-a coupling term is also due to collisional exci- 
tation. Previously, studies that have considered x-ray heating have 
ignored this effect. Recently however, Chuzhoy et at (2006) have 
pointed out that this contribution is very important and even dom- 
inates the spin-temperature value in a certain temperature range. 
In order to account for this term one should calculate the intensity 
of the Lyman-a photons due to collisional excitations, Jo. This is 
given by the following equation: 



Mr) = 



4irH(z)u a 



n H x H i{r) I a(E)J\f(E;r)dE. 



(11) 



This equation is similar to eq|7]except that instead of the fraction of 
the absorbed energy that goes to heat, /, one should use the fraction 
of the absorbed energy that goes into kinetic excitation of Lyman-a. 
The fraction, cj> a , is also parameterised by Shull & van Steenberg 
(1985) and is given by, (j> a « 0.48 (l - a; - 27 ) 1 ' 62 (where x = 
1 — xhi)- In the equation above, c is the speed of light and v a is the 
Lyman-a transition frequency. The Hubble constant as a function 
of redshift, H(z), is calculated assuming Q. m — 0.24 and Qa = 
0.76. 

The y a coupling term is (Field 1958), 



16-7T T„e /12J0 
27AioT kin m e c 



(12) 



Here, /12 = 0.416 is the oscillator strength of the Lyman a tran- 
sition, A10 is the Einstein spontaneous emission coefficient of the 
21 cm transition and e & m e are the electron charge and mass, 
respectively. 

Figure|4]shows the spin temperature of the gas for a range of 
black hole masses. The redshift and quasar lifetime (t q ) are spec- 
ified on each panel. The figure clearly shows that as the distance 
from the miniquasar increases, the temperature drops to the Tcmb 
level. The distance at which the temperature reaches the Tcmb 
asymptotic value depends on the black hole mass. For the more 
massive black holes, this distance can exceed a couple of comov- 
ing Mpc. 

The figure also shows that the maximum spin temperature is 
almost independent of the quasar mass - a detailed inspection of 
the figure shows a slight change in the maximum of T ap i n as a 
function of mass. This effect is due to the fact that the dominant 
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Figure 4. The spin temperature of the gas for a range of black hole masses. 
The redshift and quasar lifetime (t q ) are specified on each panel. The dotted 
line indicates the CMB temperature at the corresponding redshift. Note the 
different y-axis range between the 2 = 17.5 and z = 10 panels. 



coupling parameter in equation|9]around the maximum T sp i n is y a 
(by at least an order of magnitude) and is of the order of 0.01. 
Under such conditions equation [9] reduces to T sp i„ « y a Tki n , 
namely, T sp i„ tx Jo- Jo in regions where xhi <C 1 is indepen- 
dent of quasar mass as implied by the left hand side of equation [5] 
Physically, this means that when the medium is already ionized no 
additional heating of the IGM due to bound-free absorption is pos- 
sible, no matter how much radiation comes out of the quasar. 



3.2 The Brightness Temperature 

In radio astronomy, where the Rayleigh-Jeans law is usually ap- 
plicable, the radiation intensity, I(y) is expressed in terms of the 
brightness temperature, so that 



I(u) 



2iS 

„2 



hT b , 



(13) 



where v is the radiation frequency, c is the speed of light and k 
is Boltzmann's constant (Rybicki & Lightman 1979). This in turn 
can only be detected differentially as a deviation from Tcmb, the 
cosmic microwave background temperature. The predicted differ- 
ential brightness temperature deviation from the cosmic microwave 
background radiation, at the mean density, is given by (Field 1958, 
1959; Ciardi & Madau 2003), 



ST b = (20mK)(l + 5) (^-) 



Tch 



Tai 



0.0223 



10 



/0.24\ 

\n m ) 



1/2 



(14) 



where h is the Hubble constant in units of 100km s 1 Mpc 1 ,S is 
the mass density contrast, and Q m and are the mass and baryon 



© 0000 RAS, MNRAS 000. 000-000 



6 Zaroubi et al. 



30 



20 



10 



-10 
30 



20 



1 



-10 





■ - „ ; . - 




-///// \\\\\W\\ 


: t,=3 Myr VVsgvSSoSPs^^^ 


t q =20 Myr VX^sis^^^^ 


z=17.5 


z=17.5 


I H- 


+4^ H- 

- 








t q =20 Myr 


z=1 


z=10 



10 100 1000 10 100 1000 

R [comoving kpc] R [comoving kpc] 



Figure 5. The brightness temperature for the same cases shown in figure[4] 
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Figure 6. The neutral hydrogen fraction as a function of distance for 3 black 
hole masses (100, 10 4 and 10 6 MQ) for z=17.5 (left) and z=10 (right) for 
miniquasars with UV ionization energy, i.e., emitted radiation that spans the 
energy range 10.4 eV— ~ 10 4 eV. 



densities in units of the critical density. We also adopt a standard 
model universe with a flat geometry, Q,bh 2 — 0.022, Q m — 0.24 
and n A = 0.76 (Spergel et al. 2006). 

Figure[5]shows the brightness temperature for the same IMBH 
mass explored in figure[2] The curves show that the radius at which 
the differential brightness temperature is detectable increases with 
the black hole mass and the miniquasar lifetime (lefthand vs. right- 
hand panels). The maximum amplitude, however, does not depend 
on the black hole mass and depends only weakly on the miniquasar 
lifetime. This is because at the centre, T sp i„ 3> Tcmb. Hence 8Tb 
is at its maximum value which, at the mean density of the Universe, 
only depends on the redshift and cosmological parameters. 




10' 10 2 10 3 10 4 10 5 
R [comoving kpc] 




10' 10 z 10" 10* 10= 
R [comoving kpc] 



Figure 7. The differential brightness temperature for 3 miniquasars with 
black hole masses 100, 10 4 and 10 6 Mq and ionizing UV and X-ray pho- 
tons (i.e., energy range of 10.4eV < E < 10 4 eV). The quasar lifetime 
here is 3 Myr. 



4 MINIQUASARS WITH IONISING UV RADIATION 

We consider the signature of (mini-)quasars with UV radiation that 
ionizes the IGM. The different options for quasar spectral energy 
distribution has been discussed earlier. Here we follow Madau et 
al. 2004 and assume that the radiation flux spectrum is the same as 
in equation Q] except that the energy spans the range of 10.4 eV- 
100 keV. Of course in this case the quasar will ionise its immediate 
surroundings and heat up a more extended region of the IGM, a re- 
alistic spectrum will probably be between this case and the previous 
case of truncated power-law (see § refPS for a more complex en- 
ergy spectrum). Here we test three black hole masses of 100, 10 4 
and 10 6 MO at z = 10 and 17.5 with lifetimes of 3Myr. The 
10 6 Mq mass objects could be considered as progenitors of the 
SDSS 2 ss 6 quasars. The H I neutral fraction as a function of dis- 
tance from the quasar is shown in figure [6] for the three black hole 
masses at z = 17.5 (left) and z = 10 (right). 

If one assumes that the IGM is not heated relative to the CMB, 
then the quasar will heat its environment but appears as an emission 
shell around the quasar in the 21 cm brightness temperature maps. 
Figure [7] shows the differential brightness temperature around the 
same three black hole masses shown in figure [6] The clear differ- 
ence in the brightness temperature between this figure and figure [5] 
is due to the size of the ionized region around the (mini-)quasar. 



5 QUASAR FORMATION AND EVOLUTION 
5.1 Quasar evolution with redshift 

In this section we propose two very simple scenarios for the pro- 
duction and evolution of quasars at high redshift and explore the 
implications for IGM heating, ionization and the observed x-ray 
background (XRB) (Moretti et al. 2003, Soltan et al. 2003). We 
evaluate the initial mass density of black holes as a function of 
redshift, without mass accretion, with the following formation sce- 
narios: (i)- black holes as end products of stars that have formed 
through molecular hydrogen cooling, i.e., stars formed in halos 
with virial temperatures smaller than 10 4 K. (ii)- black holes 
that have been produced directly through the collapse of massive 
low angular momentum halos. In both cases, we use the Press- 
Schechter (Press & Schechter 1974) formalism with the Sheth & 
Tormen (1999) mass function to infer the number density of halos 
with a given mass as a function of redshift. 
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The mass density of black holes for the first scenario is esti- 
mated simply by calculating the number density of the most mas- 
sive halos with molecular hydrogen cooling. These are halos in the 
range of 0.1AfT 4 ^ M Mr 4 , where Mt 4 is the mass of a halo 
with virial temperature 10 4 K. This is a rough approximation for 
halos that have efficient self-shielding for H2 disassociation and can 
form pop III stars through molecular hydrogen cooling (Haiman, 
Rees & Loeb 1997a, 1997b). We henceforth refer to this scenario 
as the intermediate mass black hole (IMBH) scenario. To estimate 
the comoving mass density of the forming black holes as a function 
of redshift, we assume that at the centre of these massive halos, the 
star ends its life as a 100 M Q x {M ha i /M Ti ) black hole. The 
mass density of the forming black holes as a function of redshift 
is presented by the thick solid line shown in the upper panel of 
figured] 

For the second scenario, we estimate the number of halos 
with atomic hydrogen cooling, namely halos with virial temper- 
ature T viria i ^ 10 4 K. In order to estimate the comoving mass 
density of black holes per comoving Mpc 3 produced by this sce- 
nario, we assume that only 1% of the halos in this mass range have 
a low enough spin parameter to allow a direct collapse of the halo 
to form a massive black hole. The distribution of the spin parameter 
of halos is quite flat at the low end of the possible spin parameter 
range (Steinmetz & Bartelmann 1995), and therefore, the choice 
of 1% is rather conservative. In these halos, we take the mass that 
ends up in black holes as 10~ 3 x ^-Mhalo, where the 10 -3 re- 
flects the Magorian relation between the halo mass and black hole 
mass, and gives the baryon ratio. The comoving mass density 
of black holes produced in this type of scenario is presented by the 
solid thick line shown in the lower panel of figure [8] We refer to 
this model as the super-massive black hole (SMBH) scenario. 

To calculate the accumulated comoving black hole mass- 
density at any redshift, we assume that the black hole is accreting 
at the Eddington rate with a given radiative efficiency, e ra d- The ra- 
diative efficiency is fixed in this paper to be 10%. The cumulative 
comoving mass density is then given by the following equation, 
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p{z) = I dz p(z')e 



[M /Mpc 3 ], (15) 



where fduty is the duty cycle, which ranges from 1% to 10%, t(z) 
is the age of the universe at redshift z and tg s 0.41Gyr is the 
Eddington time-scale. 

The thin lines shown in figure [8] show the comoving black 
hole mass density as a function of redshift for several fduty val- 
ues. The calculation is done for both IMBH and SMBH scenarios. 
The case with fduty = 10% produces a black hole density rela- 
tive to the critical density of flbiackhoie(z = 6) ~ 10 -3 and 10 -4 
for the IMBH and SMBH scenarios, respectively. These values are 
too high to be compatible with the inferred black hole density at 
redshift 6. The other extreme case with fduty = 1% produces 
Quackhoie(z — 6) ~ 10 -8 for both scenarios, which is too low. 
Therefore, in the following subsections, we will focus on the results 
obtained from the cases with fduty ~ 3% and 6%. 

Recently, Begelman, Volonteri & Rees (2006) have estimated 
build up of the black hole mass density at the high redshift Universe 
that form via the 'bars within bars' mechanism. This mechanism al- 
lows for the super-massive black holes to form directly in the nuclei 
of protogalaxies, without the need for 'seed' black holes left over 
from early star formation. In their paper, Begeleman et al. (2006) 
shown black hole density as a function of redshift for two duty 
duty cycles, fduty =0.1 and 0.5. Unlike in our simple model, their 



Figure 8. Initial and evolving comoving black hole mass density as a func- 
tion of redshift. The solid thick line shows the mass density of forming 
black holes as a function of redshift. The other 4 lines show the total co- 
moving mass-density for 4 values of fduty The IMBH results are shown 
in the upper panel and those for the SMBH case are shown in the lower 
panel. 



model gives a rapid increase in the mass density of black holes until 
z ss 18 after which the black hole density evolves relatively slowly 
(roughly as logp Wacfc hole = const. + 0.086z — 0.009z 2 , which 
is obtained from a cubic spline fit to their figure 2). To summarise, 
according to Begeleman et al. the black hole density attains rela- 
tively high values early on but evolves slowly afterwords, whereas 
our model the initial density is low but the mass evolution is more 
rapid. 

Whichever the actual scenario of the evolution of black holes 
mass density in the Universe, it is clear that the mass densities ob- 
tained at high redshift contribute significantly to heating the IGM 
and decouple the spin temperature from the CMB temperature (see 
figs[5]&|7j. In the next section, we explore which of the scenarios 
we explore is consistent with the currently available observational 
constraints. 
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5.2 The Soft X-ray Background (SXRB) Constraint 

Recently Dijkstra et al. (2004) and Salvaterra et al. (2005) have 
shown that it is very unlikely that miniquasars have ionized the 
Universe without violating the observed SXRB luminosity in the 
energy range 0.5 — 2 keV (Moretti et al. 2003). In both cases 
the authors have assumed a specific black hole mass history - in- 
stantaneous in the case of Dijkstra et al. and more gradual in the 
case of Salvaterra et al. (2005). Our aim here is to check whether 
the specific black hole evolution histories proposed in the current 
study violate this observational constraint, regardless of whether 
they ionise the Universe or not. 

It will be shown that our adopted quasar duty cycle, limited 
from above by the Soltan et al. (2003) constraint, yields a diffuse 
x-ray flux that is consistent with the SXRB constraint. We assume 
a mean reionization history of the Universe according to which the 
IGM underwent a sudden reionization at redshift 6. This assump- 
tion is insensitive to our computed SXRB flux, and is conservative, 
in that it provides an upper limit on the ionizing flux from (mini- 
quasars. The SXRB is calculated for various quasar spectrum tem- 
plates. The purpose here is twofold. Firstly, to exclude from our 
models those cases that violate the SXRB constraints. Secondly, to 
explore the influence of various spectral dependences on the SXRB. 

The first template is the one we used for the quasars that have 
no UV radiation, 



where E ra n g 



10.4 eV - 100 keV. 



F(E) = A E~ a 200eV < E < lOOkeV, 



(16) 



where the calculation is made for a range of power-law indices, 
q = 2 - 0. This represents the case in which all the ionizing radi- 
ation is absorbed in the immediate vicinity of the quasar. The case 
we explored previously for the heating and ionisation fronts was 
specifically for a = 1. 

The second template, which we have also used before, repre- 
sents the case in which all the UV radiation escapes the quasar's 
immediate surroundings into the IGM. The template used here is: 



F(E) = A E' a 10.4 eV < E < 100 keV, 



(17) 



where a spans the same range as before. 

The third case we explore is the one with the template intro- 
duced by Sazonov et al. (2004) and has the form, 



F(E) 




10.4eV < E < IkeV; 
IkeV < E < lOOkeV; 
E > lOOkeV. 



Notice here that we keep the power-law index of the middle range, 
a, as the varying parameter. The reason is that quasars in the red- 
shift range 6—10 with a Sazonov et al. type spectrum contribute 
to the observed SXRB mainly in the energy range 0.5 — 2 keV. 

To proceed, we normalise the above equation with respect 
to the product of the Eddington luminosity and the radiation ef- 
ficiency, trad- This should be done at a given distance, r, from the 
quasar which we choose arbitrarily to be 1 Mpc. 

A quasar of mass M shines at e ra d times the Eddington lumi- 
nosity, namely 



L EM [M) = 1.38 x 10 38 (J^-J [erg s" 1 ]. 
Therefore A is given by: 

At*f\ trad L E dd(M) r a -1 

AIM) = ? erg s cm 

/ E~ a dE x 47rr 2 



(18) 



(19) 



In order to calculate the SXRB, we follow Dijkstra et al. 
(2004). The contribution of the soft x-ray background observed in 
the range 0.5keV < E < 2keV, given by: 



SXRB 



(l8o) 



dz cIa(z 



A{~p{z)) 



(d L (z)/Mpc)z 
T(E ' z) dE [ergs _1 cm _2 deg - ^?0) 



2(1+*) 

E~ a e 

0.5(1+2) 

In the above equation, r(E; z) represents the optical depth, 



t(E;zq) = 



dz 



(1 +z)(S/2) 



[n H i(z)a H j(-E') + n He I (z)a H e i{E')], (21) 



where E' — E(l + z)/{l + Zq), Zq is the quasar forma- 
tion redshift, tin i(z) = nn /(0) (1 + zf and n H e i(z) = 
7Wf e i(0) (1 + z) 3 are the physical density of hydrogen and he- 
lium with nn i(Q) = 1.9 x 10~ 7 cm -3 and nHe /(0) = 1.5 x 
10 -8 cm -3 . The luminosity distance, dh{z), to the black hole is 
calculated from the fitting formula given by Pen (1999) and dA is 
the angular diameter distance, 



(Ia(z) 



d L {z) 
{1 + zf 



(22) 



The division by d L accounts for the dimming of the quasar, whereas 
the multiplication by (7r/180) 2 x d\ calculates the flux received 
in a one degree 2 field of view. Moreover, the normalization factor 
is now made with respect to the mass density of black holes, and 
hence it carries an extra Mpc -3 in our units. 

Figure [9] shows the expected SXRB as a function of a for the 
IMBH and SMBH scenarios in the fduty ~ 3% and 6% cases. 
The short horizontal line at the middle of each of the panels marks 
the observational SXRB constraint. This shows that none of these 
models violate the observational constraint. The 10% case, which 
is not shown here, violates the observed constraints for almost all 
the a range. 



5.3 The number of ionizing photons per baryon 

We now calculate the number of ionizing photons per baryon emit- 
ted in the IMBH and SMBH models for the f duty = 6% and 3% 
models. The purpose of this calculation is to show that these models 
will not be able to ionize the Universe, except in the extreme case 
in which the escape fraction of the ionizing UV photons is unity 
and no recombinations take place. To estimate the number of ion- 
izing photons, one should integrate the number of emitted photons 
per unit energy over the energy spectrum of the quasars. The fac- 
tor (1 — e _T ' accounts for the absorbed fraction of photons. It also 
involves an integral over the active lifetime of the quasars down to 
redshift 6. These integrations have the following form: 

Nphotans — 47T / dz A(p(z)) fduty 



dE 



E~ a (1 - e - TW '')==- [Mpc -3 ], (23) 
E 
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Figure 9. Soft X-ray background for various spectra. The four panels show 
the SXRB level expected from the IMBH (upper panels) and SMBH (lower 
panels) scenarios with fduty = 6% (right panels) and 3%. Each panel 
shows the SXRB obtained assuming the three templates: power-law quasars 
with ionization by UV radiation (solid-lines) and without UV radiation 
(dotted-lines) and quasars with the Sazonov et al. 2004 template (dashed- 
lines). The short horizontal line in the middle of each panel marks the ob- 
servational constraint of Moretti et al. 2003. 



where dt/dz is given by, 



dt 

dz 



H (1 + z)y/(l + z) 2 (l + n m z) - z{2 + z)Q A 



[s].(24) 



Again, the mass density parameter Q m — 0.27 and the vacuum en- 
ergy density parameter Qa = 0.73. Figure[l0]shows the number of 
photons per baryon as a function of the energy spectrum power-law 
index, a. Here we note a number of features. Firstly, the maximum 
number of ionizing photons per baryon is roughly 10. This number 
is achieved in the IMBH scenario with fduty = 6% for the spectral 
templates of both Sazonov et al. (dashed line) and the power-law 
spectrum with ionizing UV radiation (solid line). Despite obtaining 
such a high number of ionizing photons per baryon, one should note 
that these two cases assume that all the quasar ionizing photons es- 
cape its immediate surroundings. Not surprisingly, the power-law 
model without ionizing photons does not produce too many ioniza- 
tions (dotted line). Note also that the number of ionizing photons 
per baryon produced by the Sazonov et al. model does not vary 
much with a. This is simply because the power-law index we vary 
in this model is in the X-ray energy range. 

Figure[TT]shows the evolution of the number of ionizing pho- 
tons per baryon with redshift. The calculation shown here assumes 
q = 1 for all three templates. The 3 left panels show results for 
the IMBH scenario, where each of the spectral scenarios is shown 
in a different panel. The right hand panels show the same for the 
SMBH case. As expected, most of the ionizing photons are pro- 
duced towards the low redshift range. The Sazonov et al. model 
produces the largest number of ionizing photons due to its steep- 
ness in the low energy range (power-law index of -1.7). Note, that 
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Figure 10. Number of ionizing photons per baryon for different spectra. The 
upper two panels show results for the IMBH scenario with the left and right 
hand panels assuming fduty = 6% and 3%, respectively. The lower two 
panels show results for the SMBH scenario with the left and right panels 
assuming fduty of 6% and 3%. The three models explored are as in the 
previous figure. 
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Figure 11. The number of ionizing photons per baryon as a function of 
redshift. The 3 left panels refer to the IMBH scenario with each of the three 
showing the number of ionizing photons per baryon for a different spectral 
template. The 3 right panels show the same for the SMBH scenario. These 
figures assume a power-law index a of 1 . 
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the fduty = 6% case produces about 10 photons per baryon at 
z = 6 normally thought to be enough to ionize the Universe and 
the same time does not violate the SXRB constraint. However, this 
model does not reproduce the Thomson r constraint (see below). 

Assuming that these curves give the actual ionization history, 
one can easily calculate the optical depth for Thomson scattering 
of CMB photons, tomb- This of course is not a self-consistent 
calculation since in order to obtain the number of ionising photons 
as a function of redshift, one has to assume an ionization history. 
This exercise is still of interest as it gives an upper limit for the 
influence of quasars on tcmb- To calculate tcmb, we assume that 
the electron density, n e = uhi x x e and x e is given by one-tenth of 
the number of ionising photons per baryon that appears in figurefTTI 
with an upper limit of unity. The tcmb found here for the IMBH 
model with fduty = 6% is about 0.03 which is only one third of 
the WMAP3 observed value (Page et al 2006). 

5.4 The predicted constraints from the Begelman et al. (2006) 
model 

As an aside, we calculate the SXRB, number of ionising photons 
per baryons and the Thomson scattering optical depth for the Begel- 
man et al. (2006) model. Here we choose the case of fduty = 0.1 
and a rigid disk model, namely, the dotted-line in the lower set of 
models in their figure 2. Since their calculation stops at z = 10 
we extrapolate their black hole mass density curves using a cubic- 
spline down to z = 6. For a Sazonov et al. (2004) type of energy 
spectrum we obtain SXRB that is about two orders of magnitude 
lower than the observed one, which is not surprising given that the 
higher the redshift of the miniquasar is, the less its soft x-ray pho- 
tons contribute to the observed background. The number of ionis- 
ing photons per baryon we obtain is about 10, usually thought to be 
the minimum number needed to ionise. However, more surprisingly 
we obtain an optical depth for Thomson scattering of, r ~ 0.075 
which is consistent with the WMAP 3 rd year results (Spergel et al. 
2006, Page et al. 2006). This result is interesting as it shows that for 
a given black hole evolution history one can satisfy all the observa- 
tional constraints and ionise the Universe solely with black holes at 
the same time. 



6 SUMMARY 

This paper explores the feasibility of heating the IGM with quasars 
without violating the current observational constraints. Such heat- 
ing is essential in order to be able to observe the 21 cm emission 
from neutral hydrogen, prior to and during the epoch of reioniza- 
tion. We have shown that miniquasars with moderate black hole 
masses can heat the surrounding IGM out to radii of a few comov- 
ing mega-parsecs. 

In this paper, two Press-Schechter based black hole mass den- 
sity evolution scenarios have been proposed, IMBH and SMBH. 
The first model assumes the black hole population is the end prod- 
uct of pop III stars that leave behind black hole masses of the or- 
der of IO-IOOA'/q. The second model assumes direct formation 
of black holes as a result of the collapse of low angular momen- 
tum primordial halos. For these two scenarios, we have explored 
three different quasar spectral templates: a power-law with ioniza- 
tion UV radiation, a power-law without ionising UV radiation and 
a Sazonov et al. (2004) type template. 

With the exception of the models that have a 10% duty cy- 
cle, we have shown that the quasars are not able to fully ionise the 



IGM - especially if one assumes the template that does not have 
ionising UV photons - while the SXRB constraint is satisfied. We 
conclude that based on the mass evolution history shown here, there 
is enough mass in the quasars to heat up the IGM by redshift 15. 
For example, for quasars with a power-law index of — 1 and no ion- 
ising UV radiation, quasars with black hole masses of 1O 3_4 M 
can heat up the IGM over a w 0.1 — IMpc comoving radius from 
the (mini-)quasar (see figure f5](. The models with 6% duty cycle 
reach such mass per comoving Mpc 3 at redshift larger than 10 for 
both scenarios. 

Curiously, for the black hole mass density evolution, with 
10% duty cycle and rigid disk model, proposed by Begelman et 
al. (2006) we find that this scenario does not violate the SXRB ob- 
servational constraint and produce about 10 ionising photons per 
baryon by z = 6, normally thought to be enough to ionise the Uni- 
verse. We also find that for a model in which the Universe ionizes 
suddenly at z=6, that this scenario predicts a Thomson scattering 
optical depth of 0.075, consistent with the WMAP 3 rd year results. 

The main result presented in this paper is "good news" for the 
new generation of low frequency radio telescopes designed to probe 
the high redshift IGM through its 21 cm emission, such as LOFAR, 
MWA and PAST. It clearly shows that the quasar population could 
easily decouple the spin-temperature from that of the CMB. 

However, since the spin temperatures achieved are not very 
high, this means that the brightness temperature will carry the sig- 
nature not only of the ionized fraction and density fluctuations, but 
also of the variations in the spin temperature. This complicates the 
interpretation of the observed brightness temperature in terms of its 
link to the cosmological fields. Nevertheless, the high spin temper- 
ature bubbles are expected to overlap before those of the ionization, 
a factor that will mitigate this complication. Furthermore, one can 
turn this around and argue that these fluctuations will teach us more 
about the ionising sources than about cosmology. An extended tail 
in the spin temperature will be a clear signature of power-law ra- 
diation, i.e., quasars, while a short tail will be a clear signature of 
thermal radiation, i.e., stars. 
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